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We study the spectrum of closed flux tubes in four dimensional SU (N) gauge theories. We do 
so by calculating the energies of the low lying states with the variational technique (whose basis 
consists of about ~ 700 operators). We study states of different values of angular momentum, 
transversal parity, longitudinal parity, and longitudinal momentum, and compare the results with 
effective string theories (ESTs) such as the Nambu-Goto (NG) model. Most of our states agree 
very well with the Nambu-Goto predictions and since most of our flux-tubes' lengths are outside 
the radius of convergence of the ESTs, then for some states it is only the NG that predicts the 
spectrum well. This strongly suggests that the ESTs can be re-summed. Nonetheless, there are 
a few states (all with negative parity and in the same representation of the lattice rotation group) 
that exhibit large deviations from the NG predictions; these deviations might provide clues to the 
nature of the effective string theory describing the large-rV QCD string. 
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1. Introduction 

We have recently shown that the closed flux tube spectrum in D = 2 + 1 SU (N) gauge 
theories can be well-approximated by the Nambu-Goto (NG) free string in flat space-time. In 
particular, such agreement was observed for flux tube lengths that are comparable to the width of 
the flux. Thus, the flux-tube approximately behaves like a fundamental string even when, naively, 
it looks more like a fat blob than a thin string. 

It is interesting to know how the closed flux tube spectrum in D = 3 + 1 SU (N) gauge theories 
behaves and how this compares to NG and other effective string theories. A pioneering attempt 
to calculate the closed flux tube spectrum in D = 3 + 1 SU (N) gauge theories has been reported 
in which investigated the spectrum for N = 3 and for a single, rather coarse, spatial lattice 
spacing of a 0.22fm. For earlier works that focused on the ground state of the closed string see 
[||] and its citations. In this paper we present results from three calculations for the closed flux tube 
spectrum. Two for N = 3 at j8 = 6.0625 (a ~ 0.09fm) and at 6.338 (a ~ 0.06fm), and one for N = 5 
at j6 = 17.630 (a ~ 0.09fm). In a further publication we will present high statistics measurements 
of the ground state for SU (3) and SU (6) at a ^ 0.09fm. All our N > 3 calculations are focused 
only on flux-tubes that carry a single unit of electric flux. 

2. Setup of the lattice calculation 

We define the SU (N) gauge theory on a four-dimensional Euclidean space-time lattice which 
is compactified along all directions with Lii x L± 1 x L± 2 x Lj sites. The length of the flux tube 
is equal to Lm, while Lj_ 1( L± 2 and Lj were chosen to be large enough so to avoid finite volume 
effects. To extract the flux tube spectrum we perform Monte-Carlo simulations using the standard 
Wilson plaquette action, S = j3 [l — ^ReTr(£/n)] . with /3 = ^?^y> and in order to keep the value 
of the lattice spacing a approximately fixed for different values of N we keep the 't Hooft coupling 
A (a) = Ng 2 (a) approximately fixed, so that j3 °< N 2 . The simulation algorithm we use combines 
standard heat-bath and over-relaxation steps in the ratio 1:4; these are implemented by updating 
SU (2) subgroups using the Cabibbo-Marinari algorithm. To measure the spectrum of energies we 
use the variational technique (e.g. see Ref. [Q] and its references). 

The closed flux tube states in D = 3 + 1 can be classified by the irreducible representations 
of the two-dimensional lattice cubic symmetry which we denote by C4. 1 It is a subgroup of 0(2) 
that corresponds to rotations by integer multiples of n/2 around the tube axis. This makes val- 
ues of angular momenta that differ by an integer multiple of four indistinguishable on the lat- 
tice, and factorizes the Hilbert space into four sectors: 7 mo d4 = 0, 7 mo d4 = dbl, / m od4 = 2. An 
additional useful quantum number is the parity P± which is associated with reflections around 
the axis ±1. Such parity flips the sign of J and so we can choose a basis in which states are 
characterised by their value of / (which can be of either sign), or by their value of |/| and P±. 
In our calculations we use the latter choice. While in the continuum states of nonzero J are 
parity degenerate, on the lattice this is exactly true only for the odd / sector. This means we 
can denote our states by the 5 irreducible representations Ai 2,£\Bi.2 of C4 whose J and P± as- 
signments are: {A\ : |7 mo d4| = 0,Pj_ = +} ,{A 2 : |/ mo d4| = 0,Pj_ = -}, {E : \J mo d4\ = h p ± = ±}, 

'Since the global rotational symmetries of a flux-tube wound around a compact direction are those of a two- 
dimensional space. 



2 



Spectrum of closed fundamental flux tubes in D = 3 + 1 



Andreas Athenodorou 



{Si : |7 mo d4| = 2,.Pj_ = +}, and {B2 : |7 mo d4| = 2,Pj_ = — }. All the representations of C4 are one- 
dimensional except for E which is two-dimensional. 

Two additional useful quantum numbers include the longitudinal momentum pii carried by the 
flux-tube along its axis (which is quantized in the form p\\ = 2nq/L^;q G Z) and the parity Pu with 
respect to reflections across the string midpoint. Since Pu and pu do not commute, we can use both 
to simultaneously characterise a state only when q = 0. Also, since the energy does not depend on 
the sign of q, we only focused on those with q > 0. 

The operators we construct have shapes that lead to certain values of J,P±,P^, and q. This 
is achieved by choosing a linear combination of Polyakov loops whose paths consist of various 
transverse deformations and various smearing and blocking levels (again see [Q]). All the paths 
used for the construction of the operators are presented in Table [j] and all together form a basis of 
around 700 operators. Let us show how to construct an operator with a certain value of / mo d4: begin 
with the operator a that has a deformation extending in angle a within the plane of transverse 
directions. We can construct an operator 0(7) that belongs to a specific representation of C4 by 
using the formula: 0(7) = En=i,2,3,4 ei/ni 0n§ • It 1S straight-forward to show that 0(0) belongs to 
either A 1 or A2 (depending on its value of P±), that 0(1) belongs to E, and that 0(2) belongs to B\ 



or #2- The projection onto certain values of P± and Pu is demonstrated pictorially in Eq. ( |2.1[ ) for 
an operator of 7 mo d4 = 0. 




=Tr 



If i = j = k = +1 then the operator projects onto {A 1 , 7^ = +}, if / = + l,j 
projects onto {A2,P\\ = +}, if i = — l,j = +l,k = — 1 it projects onto {Ai,Pii = 
i = j = — l,k = +1, it projects onto {A2,P\\ = — }. 



(2.1) 
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Table 1: All the transverse deformations used for the construction of the operators. 

3. Theoretical expectations for the spectrum from effective string theories 

Let us first think about the flux-tube as a string of length / = aLu winding around the torus. 
The classical configuration of the string spontaneously breaks translation symmetry and so we 
expect a set of Nambu-Goldstone massless bosons to appear at low energies. These bosons reflect 
the transverse fluctuations of the flux-tube around its classical configuration, and in D space-time 
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dimensions there are D — 2 of them. Below, we describe the current theoretical predictions for the 
excitation spectrum of these bosons. 

3.0.1 The Nambu-Goto model 

The Nambu-Goto model describes relativistic strings [ph. The action of the model is propor- 
tional to the area of the world sheet swept by the string as it propagates in time. This model is 
self-consistent quantum mechanically only in D = 26 dimensions (see for example second refer- 
ence in []|]), but there are claims in the literature that, for any value of D, this model can serve as 
an effective low energy field theory for long strings From here on we refer to such low energy 
theories as Effective String Theories (ESTs). The spectrum of the Nambu-Goto model for D = 4 is 
given below (e.g. see second reference in Ref. [§]). 



E NG (l) = ^(ol) 2 +47io(N++N L +n++Nz-\) + (^) ■ (3.1) 

Here a is the string tension and and Ng are the occupation numbers of the bosons that move to 
the left and to the right (the ± superscript denotes the spin that they carry). This means that each of 
these occupation numbers is defined to count the energy units carried by the bosons: N = YX=i knk 
(here is the number of bosons carrying momentum k). The net longitudinal momentum carried 
by the bosons is given by q = N$ + N£ — N£ — , and the net angular momentum of a state is 
given by J = +Nr — N£ —Ng. It is useful to make a connection to Regge theory by writing 
Y,i=±lLj=L,RNj = J + m and then interpreting the integer m as counting the daughter trajectories of 
a certain angular momentum J. As usual, leading Regge trajectories of angular momentum / can 
be degenerate in energy with daughter trajectories of lower spin states. Below we refer to states 
that have the same value of J + m as 'being in the same NG level'. 

Since we think of the NG model as an EST, which may be justified only for long strings [0], 



we can expand Eq. (3.1) for l^Jo 3> 1 and we get (for simplicity we set q = 0, and follow the 



convention in denoting / + m by n/2) 

An f D-2\ 8tt 2 / D-2\ 2 32n 3 ( D-2\ 3 . 7n 

£N0(/) = CTi+ _(„__j __(„__) +_(„__) +0(r >). (3.2) 

We note in passing that because the NG model is only one possible candidate of an EST, then in 
the language of effective field theories, it may differ from other candidate ESTs by the values of 
certain low energy constants (LECs). These LECs would make the most general EST spectrum 



differ from Eq. ( p^2| ) by the coefficients in the \/l expansion. In the next subsection we discuss a 
systematically controlled approach to construct the most general EST of the QCD flux-tube. 

3.0.2 Effective string theory approaches 

A systematic EST study that would describe the QCD flux-tube was pioneered by Liischer, 
Symanzik, and Weisz in Ref. Such an EST approach produces predictions for the energy of 
states as an expansion in 1//. Terms in this expansion that are of 0(1 /l p ) are generated by (p+ 1)- 
derivative terms in the EST action whose coefficients are a priori arbitrary LECs. Interestingly, 
these LECs were shown to obey strong constraints that reflect a non-linear realization of Lorentz 
symmetry ^ ^, [K)|] , and so to give parameter free predictions for certain terms in the 1 // expansion. 
We review these predictions below. 
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First, since we focus on closed strings, p can only be odd (terms in the energy that come with 
even powers of 1 // appear only if there are boundary terms in the action of the EST and so do 
not exist for closed strings). The analysis with p = 1 was performed in Ref. [{7|] and was shown 
to yield the 0(1/1) 'Luscher term' in the ground state energy, whose universal coefficient depends 
only on D. Since the 2— derivative action is a free theory, then its spectrum of excited states is 
that of (D — 2) massless bosons (the Luscher term is the zero point energy of these bosons). The 
predictions of Ref. [Q] were that to 0(1/1) the flux-tube energy is given by the first two terms in 
right hand-side of Eq. ( |3.2[ ). 

The 4— derivative terms were analysed in Ref. []|] and for D = 3 were shown to yield an 
0(l/P) term that is identical to the third term in Eq. (3.2). Ref. [1C] showed that this matching 
between the 0(1 /I 3 ) term in the NG prediction and in any general EST holds also for D = 4. 
Ref. [[K]] also analyzed the 6— derivative terms and showed that for D = 3 they yield the fourth 
term in the r.h.s. of Eq. (3.2), while for general states in D = 4, the coefficient of the 0(1 /I 5 ) term 
may differ from the one in Eq. (3.2). Nonetheless, the energy of n = state in the D = 4 case 
is special and Ref. [1C] showed that its 0(l/l 5 ) term is indeed given by Eq. (3^), and that the 
0(1/ 1 5 ) term in the average over the energies of states that are in the same NG level is identical to 
the 0(1 /I 5 ) in Eq. Q. 

A different approach to EST was proposed by Polchinski and Strominger in Ref. [11]. Tech- 
nically it uses a different gauge fixing of the embedding coordinates on the world-sheet (conformal 
gauge instead of the static gauge choice used in the ESTs following Ref. [0]). Here the constraints 
obeyed by the LECs allows one to maintain the conformal symmetry of the world-sheet even out- 
side the critical dimension of D = 26. Ref. [11] showed that as a result of these constraints, the 
0(1/1) term in the 1 // expansion is the same as in the NG model — it is given by the Luscher term 
appearing in Eq. ( |3.2[ ). Much more recently, Drummond [|l^] showed that the 0(1 /I 3 ) term is also 
identical to the one appearing in the l/l expansion of the NG model (Eq. (3.2)). Finally, the recent 
Refs.[ 13] claim that even higher order terms in the l/l expansion are identical to the corresponding 
terms in Eq. (3.2); these claims seem to contradict the results of Ref. JIo| ] for D > 4. 

4. Results 

4.1 Ground states in the q = 0, 1, 2 channels 

In Figure [j] we present a comparison of the q = 0, 1 , 2 NG ground states with our data: SU (3) 
and a ~ 0.09fm in the left most plot, SU (3) and a ~ 0.06fm in the central plot, and SU (5) and 
a ~ 0.09fm in the right most plot. The NG ground state for q = is the state of no bosons, and 
has trivial quantum numbers {J = 0,P± = +,P\\ = +}■ Indeed our data shows that the ground state 
is in the A\ representation with q = and Pu = +. We use the measured energy of this state to 



extract the string tension by fitting it with the ansatz E gs (l) = ^j(ol) — ^y- + ; we find 

C ~ 0(1 — 10) and that for l^fa > 3 the correction term is negligible. 

The NG ground state in the q = 1 channel has 7=1 and is P± degenerate, and so we measure 
only the energy of the P± = + channel in the E representation. This measurement is presented in 
Fig. [j] where we see that the NG prediction is in agreement with our data. Finally, NG predicts that 
the ground state for q = 2 should be five-fold degenerate, consisting of a state with {J = 0,P± = +}, 
two states with {J = l,P± = ±}, and two states with {J = 2,P± = ±}. We find this to be consistent 
with the ground states in the A\ ,E and B\ 2 representation (again, as for q = 1, we measure only the 
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P± = + states in the E representation). Comparing the three plots in Fig. [j] we see that the 0(a 2 ) 
and 0(1 /N 2 ) corrections of our data are small compared to our statistical errors. Thus, within our 
level of accuracy, the agreement of our data with the NG model is largely insensitive to the lattice 
spacing and \/N 2 corrections. 
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Figure 1: Energies of the lightest states that correspond to q 
4.2 Further q = states 

We now turn to examine the first excited state in the q = sector and present the results in 
Figure |2[ In the NG model, this energy level is four-fold degenerate with a single state from each 
of the following representations: {J = 0,P ± = P\\ = +}, {J = 0,P ± = P\\ = -}, {J = 2,P ± = P\\ = 
+}, and {J = 2,P± = — ,P\\ = +}. On the lattice, this would imply that states from (A\,P^ = 
+), (A2,P\\ = —),(Bi,Pi\ = +), and (B2,P\\ = +) are all degenerate (up to 0(a 2 ) corrections). 

Instead, what we find is that while the states belonging to {A\,P\\ = +),(B\,P\\ = +), and 
{B2,P\\ = +) are all quite close to each other and to the NG model, the state in (A2,P\\ = — ) 
is anomalously different and shows substantial deviation from NG. It is tempting to expect that 
this state's energy would eventually approach the NG prediction (and is perhaps reflecting a large 
coefficient multiplying the term that controls this deviation from NG in the EST). Nonetheless, 
an equally likely possibility is that it does not converge to NG. In fact, observe that the energy of 
this state is higher than the ground state energy by approximately an equal amount throughout the 
distance range that our simulations are able to explore. Therefore, at large enough /, it might cross 
the NG prediction, as a massive state would. We do not know which possibility provides a better 
explanation for the strikingly 'anomalous' way that the energy of this state behaves. This behaviour 
does not change as we decrease the lattice spacing (central plot) or increase (right-most plot). 

In Figure § we also provide a comparison of the four states with expansions of the Nambu- 
Goto square root order by order in terms of 1/7 up to 0(1 /I 5 ) - see Eq. (^2) (note that according 
to Ref. [10] it is only up to 0(l/P) that we can trust Eq. ( |3.2[ )). Excluding the anomalously 
behaving ground state in the {A2,P\\ = — } channel, the other three states are obviously better 
described by NG than by any other EST. This reflects a simple fact: nearly all our data is beyond 
the radius of convergence of the 1 // expansion (which can be estimated from the EST series to be 
(/- v /a) converge ~ 3.5 for this particular state). Thus, the fact the our data is so close to NG even 
beyond the radius of convergence suggests that the EST approach can be somehow resummed to 
yield the NG square root of Eq. (3T) plus some small corrections. In our 2 + 1 study [jl|] we also 
saw an agreement with NG beyond the radius of convergence of the EST. 
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Figure 2: Energies of the four states with q = expected to be characterised by n = 1 . 



4.3 Further q = 1 states 

We now discuss further results for states with q = 1 — see Figure [| Here the first excited 
NG energy level (n = 3/2) should be ten-fold degenerate. These states are the ground states of 
{J = 0,P ± = +}, {J = 0,P ± = -}, {J = 3,P X = ±}, {J = 2,P X = +}, {J = 2,P ± = -} and the 
first excited states of {/ = 1,P± = ±}. On the lattice these states fall into two degenerate pairs 
of states in E, and four more states that belong to A\£ and Bi 2- The parity degeneracy in the E 
representation allows us to calculate only the P± = + states and we therefore expect six degenerate 
states in the NG model. 

Unfortunately, our basis of operators was insufficient to successfully isolate all these six states, 
and we were able to extract only four of them. We found that the energy of the q = 1 lowest energy 
states in the A\ and B\ representation, and second lowest energy states in the E representation, 
agree fairly well with the energy of the n = 3/2 NG level. In contrast, the ground state of q = 1 
in the A2 representation (which is naively associated with the / = 0, P± = — channel) appears 
anomalous: it has a large deviation from the NG curve and does not show any sign of convergence. 
This is true also on our finer lattice (central plot) and for SU (5) (right-most plot). 
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Figure 3: Energies of the lightest five distinguishable states with q = 1, 
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5. Summary 

We calculated energies of 13 states in the spectrum of closed strings in 3 + 1 dimensional 
SU (N) gauge theories. We study the gauge groups SU (3) (with lattice spacing a fa 0.09,0.06fm) 
and SU (5) (with a as 0.09fm) and the string lengths in the range 0.9fm < I < 1.8fm. 

Most of the energies we measured show a convincing agreement with the Nambu-Goto (NG) 
model. We compare our results with recent predictions from other effective string theories (ESTs), 
but because most of the flux-tubes that we probe have lengths that are beyond the radius of conver- 
gence of the ESTs (for the excited states), this comparison fails. The only model that describes our 
data well, even for short flux-tubes and excited states, is the NG model. This was seen also in our 
previous work on 2 + 1 dimensions (see Ref. [jl|]), and seems to suggest that the EST expansion can 
be resummed. Despite the overall good agreement with NG, we do find large deviations for certain 
states that have negative parity and that belong to the A\ representation of the lattice group (which 
naively corresponds to zero angular momentum). We see these deviations also for states that carry 
both zero and one unit of longitudinal momentum. 

There are many avenues one could take on the lattice to make progress towards establishing 
what is the effective string theory of the QCD flux-tube. These future studies may include the search 
for massive excitations like breathing modes, attempting to accurately test the current theoretical 
predictions (see Section ||) within their radius of convergence, studying the open string spectrum, 
etc. We also look forward to theoretical progress that would allow one to understand the results we 
presented in this proceedings. For example, how can the flux-tube behave like a NG string below 
the radius of convergence of the effective string theory expansion? and what makes the states in 
the Ai representation have large deviations from the NG model? 
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